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^ Abstract 
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A method of Feynman diagrams summation, based on using Schwinger-Dyson 
equations and Ward identities, is verified by calculating some four-loop diagrams 
in N = 1 supersymmetric electrodynamics, regularized by higher derivatives. In 
particular, for the considered diagrams correctness of an additional identity for 
Green functions, which is not reduced to the gauge Ward identity, is proved. 



1 Introduction. 



Current indirect proofs of existence of the supersymmetry in the Standard model 
make the problem of calculation of quantum corrections in supersymmetric theories 
especially urgent. Due to the supersymmetry the ultraviolet behavior of field theory 
models is essentially improved. For example even in theories with unextended super- 
symmetry it is possible to suggest the form of the /3-function exactly to all orders of the 
perturbation theory. In the case of N = 1 supersymmetric electrodynamics, which will 
be considered in this paper, this /3-function (that is called the exact Novikov, Shifman, 
Vainshtein and Zakharov (NSVZ) (5- function) is pQ: 

/8(a) = ^(l- 7(a)), (1) 

where 7(a) is the anomalous dimension of the matter superfield. 

Explicit calculations, made with the dimensional reduction confirm this pro- 
posal, but require a special choice of the subtraction scheme [Hill]- Explicit calcula- 
tions in two- Ej and three-loop [Zj approximations for the N — 1 supersymmetric 
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electrodynamics with the higher derivative regularization jHJ |U] reveal that renormal- 
ization of the operator W a C ab Wb is exhausted at the one-loop and the Gell-Mann-Low 
function coincides with the exact NSVZ /3-function and has corrections in all orders of 
the perturbation theory. 

An attempt to perform explicit calculations exactly to all orders of the perturbation 
theory was made in Ref. [TU]. According to this paper it is possible to calculate a large 
number of Feynman diagrams exactly to all orders of the perturbation theory using 
Schwinger-Dyson equations and Ward identities. Nevertheless, some diagrams can 
not be obtained by this method. However, explicit calculations (up to the three-loop 
approximation inclusive) show, that the sum of undefined contributions is 0. This 
means existence of an identity, which is not reduced to the Ward identities and can 
be graphically written in the form (In the analitic form this identity is given 

by Eq. (JHUjO . It relates the two-point Green function of a matter superfield and a 
three-point Green function, in which external lines correspond to a gauge field, a chiral 
matter superfield and an external source introduced in a special way. This source is 
an arbitrary scalar superfield. 

Identity (JHU|) is important because it restricts nontrivially the structure of the di- 
vergences of the theory. Really, due to the supersymmetry and gauge invariance the 
two-point Green function of the gauge field should be proportional to Vd 2 Ui/ 2 V, the 
coefficient being an unknown function. (It is a supersymmetric analog of the transver- 
sality condition in the ordinary electrodynamics.) Using the Schwinger-Dyson equa- 
tions and Ward identities it is possible to rewrite this function via the two-point Green 
function of the matter superfields and another function, which is not fixed by the gauge 
invariance. Actually this means, that the two-point Green function of the gauge super- 
field is equivalently rewritten in a different form, as earlier, up to an undefined (from 
the Ward identities or, equivalently, from the gauge invariance) function. However, 
according to new identity the terms, which depend on the unknown function, are ac- 
tually equal to zero. So this identity really removes the arbitrariness in the structure 
of the divergence. 

Certainly, it is highly desirable to find a true reason of this identity - for example, 
some symmetry. But first it is necessary to verify, that it really takes place and is not 
caused by an accidental coincidence. 

The considered identity is nontrivial starting from the three-loop approximation. 
In Ref. ^T] it was proved for a special class of diagrams exactly to all orders of the 
perturbation theory. Such diagrams contain the only loop of the matter superfields 
and any two cuts of this loop do not make the diagram disconnected. (In this case the 
technical part of the proof is simpler.) Possibly in the general case the proof can be 
made similarly. In order to confirm that such a proof is possible, it is desirable to verify 
identity ()3U|) for other diagrams. This requires making explicit four-loop calculations. 
These calculations are made in this paper for a special group of four-loop diagrams, 
which can be made disconnected by two cuts of the matter superfields loop. 

The paper is organized as follows: 

In Sec. |2] the basic information about N = 1 supersymmetric electrodynamics 
and its regularization by higher derivatives is reminded. The method of summation 
Feynman diagrams, based on using Schwinger-Dyson equations and Ward identities is 
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described in Sec. 01 In this section we write an identity, to which the Green functions 
are supposed to satisfy. Four-loop verification of the results, presented in Sec. El is 
made in Sec. EJ A brief summary of the paper is presented in the Conclusion. 



2 N = 1 supersymmetric electrodynamics and 
higher derivative regularization. 



The massless N = 1 supersymmetric electrodynamics with the higher derivatives 
term in the superspace is described by the following action: 



4e 2 



d*xc?6W a C ab (l 



d 



2n 



A 2n 



W b + - d 4 xd 4 9 



L*„2V . 



+ 0*e- 2V > . (2) 



Here and are the chiral matter superfields and V is a real scalar superfield, which 
contains the gauge field as a component. The superfield W a is a supersymmetric 
analog of the stress tensor of the gauge field. In the Abelian case it is defined by 



W a = -D(l- l5 )D (l + l5 )D a V 
lb L 



(3) 



where D is a supersymmetric covariant derivative. It is important to note, that in the 
Abelian case the superfield W a is gauge invariant, so that action (j2J) will be also gauge 
invariant. 

Quantization of model (j2J) can be made by the standard way. For this purpose it is 
convenient to use the supergraphs technique, described in book [12J in details, and to 
fix the gauge invariance by adding the following terms: 



s " MB (^i 1 +d ^) v+ vd2d \ i +d Pi v \ « 

where 

D 2 = ^D(l + l5 )D; D 2 = ^D(l- l5 )D. (5) 

After adding such terms a part of the action, quadratic in the superfield V, will have 
the simplest form 

s„ + s gf = —J d'x d'e vd 2 [i + —)v. (6) 

In the Abelian case, considered here, diagrams containing ghost loops are absent. 

It is well known, that adding of the higher derivative term does not remove diver- 
gences in one-loop diagrams. In order to regularize them, it is necessary to insert in 
the generating functional the Pauli-Villars determinants |13j . 
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The generating functional can be written in the form 

Z= f DVD ( j)D&]j(detPV(V,M i )y\xp(i(S ren + S g f + S s + S l p )). (7) 

Here 

S ren = ~^Z 3 (e, A/fj,) Re J d 4 x d 2 6 W a C ab (l + ^) W b + 

+Z(e, A/ n) i J d 4 x d 4 6 (fe 2V (j) + 0*e" 2V 0) (8) 

is the renormalized action, the action for the gauge fixing terms are given by Eq. (j^J) (It 
is convenient to substitute e by the bare charge e in it), the Pauli-Villars determinants 
are defined by 

(detPV(V,M)) 1= / exp (iSw) , (9) 

where 

5 Py = Z(e,A/fi)^ J d 4 xd 4 6 ($*e 2V $ + <J* e - 2y l>) + 

+~ y d 4 xd 2 6M$$ + ^ J d 4 xd 2 9M®*<S>\ (10) 

and the coefficients q satisfy conditions 

$> = 1; ^ Ci M 2 = 0. (11) 

Below we will assume, that Mj = OjA, where are some constants. Insertion of Pauli- 
Villars determinants allows to cancel remaining divergences in all one-loop diagrams, 
including diagrams, containing insertions of counterterms. 
The terms with sources are written in the form 

S s = Jd 4 xd 4 0JV + J d 4 xd 2 6(j<P+]<P) + J d 4 xd 2 9 (j*<f>* +]*4>*y (12) 
Moreover, in generating functional (J7J) we introduced the expression 

S 4>o = \J d 4 x d 4 6 (0* e 2V <P + 0* e 2V 4> + 0* e" 2y + 0* e~ 2V o ) , (13) 

where 0o, 0o, 0o and 0q are scalar superfields. They are some parameters, which are 
not chiral or antichiral. In principle, it is not necessary to introduce the term S$ into 
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the generating functional, but the presence of the parameters 0o is highly desirable for 
the investigation of Schwinger- Dyson equations. 

In our notations the generating functional for the connected Green functions is 
written as 

W = -i]nZ, (14) 
and an effective action is obtained by making a Legendre transformation: 

r = W - J d 4 xd 4 9JV - J d 4 x d 2 9 (j <f> + j 0) - J d 4 x d 2 6 + , (15) 

where the sources J, j and j is to be eliminated in terms of the fields V, 4> and (ft, 
through solving equations 

v= jj ] 0= ir 0= tt ( 6) 

3 Summation of Feynman diagrams in TV = 1 super- 
symmetric quantum electrodynamics 



From generating functional (J7J) it is possible to obtain the Schwinger-Dyson 
equations, which are graphically written as 




^ — 'WWJ^ (lAAAA, -\- 

^a^ 1 la M A/X/ (17) 

(2) 

where T v is a two-point Green function of the gauge field. We will always set the 
renormalization constant Z — 1, because the dependence on Z can be easily restored 
in the final result and this dependence is not essential in this paper. 
The double line in Eq. (fT7|) denotes the exact propagator (Z = 1) 



^ ' x Sl y , (18) 



Scj)* x 5(j) y j Ad 2 G 
in which the function G(q 2 ) is defined by the two-point Green function as follows 

5 2 T DID 2 



7^(^)5" (19) 



where 5® y = 5 4 (x — y)5 4 {6 x — 9 y ), and the lower indexes denote points, in which con- 
sidered expressions are taken. 

The large circle denotes the effective vertex, which is written as 
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= d 2 u 1/2x [Dy xy Dy xz )F( q ") + 

p=0 

+^G\q 2 )D^ 5 D x (Dl5l y Dl5l z ) +\Dl8l y Dl% z G{<l% (20) 



due to the Ward identities. The strokes denote derivatives with respect to q 2 



n V2 



1 



D a D 2 D n 



16d 2 



1 



— —-D a D 2 D a 
16d 2 



(21) 



is a supersymmetric transversal projector, and F(q 2 ) is a function, which can not be 
defined from the Ward identities. Here 



D a = 
D a = 



2 (! +^)D 



(22) 



Two adjacent circles denote an effective vertex, consisting of 1PI diagrams, in which 
one of the external lines is attached to the very left edge. Such vertexes are given by 

DDI 



6 2 T 



1 6 I 2 5 



—cm* 



in the case of one external VMine (the vertex in the first diagram of Eq. ()17Jl ) and 



(23) 



6 3 T 



6V 



X U^yU^ Qz 



p=0 



1 6 6 

4 6V X 50, 



exp 



2 _6_ 

i 5J 2 



+ 2V 



-2d 2 U 1/2x ( D 2 Jl y 6» xz )F(q 2 ) + -D«C ab D 2 x ( D^D'X, ) f(q 2 ) + 



p=0 
S2r8 n&x8 



'xy 



-^G'(q 2 )Dri 5 D x (DX v 6U -^DXAG(q% (24) 
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1 = 



'xy- 



in the case of two external ^-lines (in the second diagram of Eq. (|17)l). Here f{q 2 ) is 
one more function, which can not be found from the Ward identity. 
The expression 



d r (2) 
dlnA v 



(25) 



p=0 



can be calculated jTU] by substituting solutions of Ward identities into the Schwinger- 
Dyson equations. 
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A result of calculation of any diagram can be written in the form 

/ d * 6 { Vd2U ^ VA (p 2 ) + VBtfj), (26) 

where A and B are some functions. It is actually (with a small modification) a su- 
persymmetric analog of expansion to the longitudinal and transversal parts. Due to 
the Ward identity the sum of diagrams contains only transversal parts, i.e. terms, 
proportional to Vd 2 Ui/2V . It is confirmed by calculations, made in Ref. ^U] : Terms, 
proportional to V 2 , are in the sum of diagrams, and we will not write them here. 
Terms, proportional to Vd 2 Ui/2V , in the effective diagrams are given by 

— ^/^(_-_ v). (28) 

(Here we omit integration with respect to the external momentum and external 9.) 
Expressions for these diagrams also contain terms, proportional to V 2 and similar 
terms with the Pauli-Villars fields, which we omitted for brevity. Their form can be 
found in Ref. jTU]- Terms with the Pauli-Villars fields give additional contributions, 
which cancel the remaining divergences. Investigation of them is a bit more complicated 
and requires making similar calculations in the massive theory. We will always assume 
existence of terms with the Pauli-Villars fields, but will not write them explicitly. 

Note now, that terms, containing the unknown function F, are completely cancelled 
in the sum of contributions (}2~Tj) and ([28)1 . However there are terms, containing the 
unknown function /: 



d 6 2 AT 



din A 5V x 5V y 



p=0 



"K^Jw? (w4 Hq2G2) -^- (pv) \ (29) 



where (PV) denotes similar terms with the Pauli-Villars fields. Nevertheless, explicit 
three-loop calculations show that the following identity takes place: 

d r *< /( « 2 > -o, (so) 



din A 7 (2tt) 4 q 2 G(q 2 ) 

where the function G is defined from the two-point Green function of the matter 
superfield by Eq. ([19)1 . and the function / - from the three-point Green function with 
the external lines V, 4> and 0q by Eq. (JUJ). 

An attempt to prove identity ()30|) exactly to all orders was made in Ref. jTTJ by 
approximately the same method, which is used for the proof of the Ward identity 
diagram by diagram. Although application of this method appears to be possible, the 
proof was not made in the general case due to the technical difficulties. It was made 
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only for diagrams, which contain the only loop of the matter superfields and remain 
connected after any two cuts of this loop. In the general case the proof, made by the 
same method, seems to be possible, but much more complicated technically. 
Identity (J3U|) can be graphically written as 




0, 

(31) 

D a V X Vy 

that was also proved in Ref. [TTj. The symbol (Ii) a (here we use notations of Ref. [H] ) 
means, that the double line in this case corresponds to the expression 

' D°D 2 JL (32) 



2d 2 G(d 2 j 

instead of the effective propagator (fT8|) . which is proportional to D 2 x D 2 x b\ z . Note, that 
equality (|3T)|) is not valid in the massive case. Its corresponding modification can be 
found in [TU] . 



4 Four-loop calculations 



Note now, that expressions ()27|) and (|28|) allow to find not only sums of all diagrams 
with two external lines of the gauge field, but also sums of special classes of such 
diagrams. Such classes of diagrams are obtained from a frame, to which external lines 
are attached by all possible ways. For example, in Fig. Q we present a frame for 
diagrams, which are investigated in this paper. The dots denote all possible points, to 
which two external lines of the gauge field can be attached. 




Figure 1: Schematic picture of the class of diagrams, considered in the paper 

If a frame contains the only loop of the matter superfield and can not be made dis- 
connected by two cuts of this loop, then for the corresponding class of diagrams identity 
(131 \ was proved in Ref. [H]. However it is desirable to verify, if this identity takes 
place in the other cases. For this purpose we consider diagrams, which are obtained 
from the frame, presented in Fig. Q by adding two external lines. (Such diagrams can 
be evidently made disconnected by two cuts of the loop of matter superfields.) 

We will calculate the considered diagrams by two ways: 

1. using Eqs. (|2"Tj) and (j2Hl) and the functions G, f and F, obtained preliminary. 
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2. by explicit calculation using the supergraph method. 

Thus it is possible to perform a four- loop verification of Eqs. (12 7|) and (|28jh and 
also identity (|30|). 

In order to find diagrams which will be essential for finding the unknown functions 
G, f and F in the considered case, it is convenient to use the following simple specu- 
lations: A diagram, presented in Fig. Q can be considered as a formal product of one- 
and two-loop diagrams with <fi* and <ft external lines (pairs of their ends are identified): 



/0\ 




or as a three-loop diagram with the identified ends: 




The parts of the diagram, obtained by this way, are the Feynman diagrams for finding 
the function G. To define the functions / and F it is necessary to add one more 
external line of the superfield V to such diagrams. 




Figure 2: The diagrams, used for calculation of the function G 



Thus it order to find the function G in the considered case it is necessary to calculate 
diagrams presented in Fig. |21 Then the function G is obtained according to definition 
(Unj). One-, two- and three-loop parts of the function G, defined by diagrams in Fig. 
I2J will be denoted by G\, G2 and G 3 respectively. Expressions for them are presented 
in Appendix [X] The complete function G (certainly without diagrams, which are not 
essential for this paper) in the considered approximation is given by 

G(q 2 ) = 1 + G 1 + G 2 + G 3 , (33) 

where the unity is a tree contribution. 

The functions F(q 2 ) and f{q 2 ) are defined from three-loop Green functions by 
Eq. (124)1 . As in the case of the function G, in order to obtain the contribution, which 
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corresponds to the considered class of the four-loop diagrams, it is sufficient to calculate 
only diagrams of a special form. These diagrams are obtained from diagrams, presented 
in Fig. |21 by all possible insertions of one external ^/-line. 

We will denote one-, two- and three-loop contributions to the function / by fi, fa 
and /3 respectively. Similar notation we will use for the function F. Because in the 
tree approximation these functions are and fa = 0, in the considered order we have 

F = Fi + F 2 + F 3 ; / = fa + fa. (34) 

Expressions for fa, fa, fa, Fi and F 2 , obtained by calculation of the above pointed 
diagrams, are presented in Appendix |XJ An expression for F 3 has been calculated, but 
it is not presented because it is very large. (It is not required for verification of identity 

Using the obtained expressions it is possible to verify identity (J3UJ). With the 
considered accuracy we have 




Substituting here expressions for the functions G\, fa and fa from Eqs. (jHSJ), dlSj) , Q44j) 
we obtain, that the integrand can be written as a total derivative with respect to the 
momentum q: 



d*q f 



(2tt) 4 q 2 G 
l + (-l) n 



le 



d 4 qd*kdHd A r 



(2 



7T 



16 



„2n 



i + (-ir 



,2n 



A 2n 







{2q + k + iy 



i + (-ir 



l2n 



A 2r 



dq n y k 2p r 2 q 2( k + q y( q + iy( q + r y( k + q + iy 



(36) 



This equality can be verified by calculating the derivative with respect to q^ using the 
Leibnitz rule and comparing the result with Eq. (J35|) . in which expressions for the 
functions G\, fa and fa are obtained by explicit calculation of diagrams. 

Because the integrand in Eq. (}3l)j) is a total derivative with respect to g M of the 
expression, which goes to in the limit q — > oo, expression (}36|) is 0. This means that 
identity (jHOj) is correct in the considered approximation and for the considered class of 
diagrams. 

Note, that in this case using of the higher derivative regularization is essential. In 
spite of the parameter A is encountered only in propagators of the superfield V, which 
do not contain the momentum q, this regularization allows to make calculations in the 
limit p — > 0. It follows from existence of limit (J25j) . which was proved in Ref. |lUj . 
Taking this limit is senseless without the higher derivative regularization, because this 
produces infinities. Another important reason of using the higher derivative regulariza- 
tion is the possible problems with integration with respect to the momentum q in Eq. 
(J29|) in D ^ 4 dimensions. For example, it is possible, that in D dimensions expression 
(p3H|) can not be presented as a total derivative. In principle, identity (f3T)|) can be con- 
sidered as a formal relation in the non-regularized theory, which is made sensible by 
the regularization. However, most likely it is not valid for an arbitrary regularization in 
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the regularized theory. Because supersymmetric theories are usually regularized either 
by the dimensional reduction or by the higher derivative regularization, it is sensible 
to say, that identity (j3*Uj) is valid in the theory, regularized by higher derivatives. 

It is also desirable to verify the method of summation of Feynman diagrams by using 
the Schwinger-Dyson equations and Ward identities. For this purpose it is possible to 
calculate both effective diagrams in Eq. (J 17)) in the four-loop approximation explicitly 
and compare the result with Eqs. (|27j) and (|28|) . In considered approximation 

7^ ln < G2 )^^( Gs - G '' G2 ) ; < 37 > 

F If \ 
-a^ = -2\F 3 - F 2 ■ Gi - F x ■ G 2 ). (38) 

Explicit calculations were made by the method, proposed in Ref . |T3] , which simplifies 
finding a part of diagram, proportional to Vd 2 Ui/2V . Nevertheless, it does not allow 
to obtain a part, proportional to V 2 . That is why the verification of Eqs. ()27|) and (|28|) 
was made only for a part, proportional to Vd 2 Ui/2V . In the both cases this verification 
completely confirms them. 

As a small technical remark let us note, that making this verification it was nec- 
essary to take into account, that a large number of ordinary Feynman diagrams con- 
tributed both to the first effective diagram in Eq. ()17j) and to the second one. For 
example, it is easy to see, that a contribution of diagram, presented in Fig. |3] 




Figure 3: One of the diagrams, giving contribution to both effective diagrams 

is divided into parts 3/4 and 1/4, which correspond to the first and to the second 
diagrams in Eq. (JT7J). 



5 Conclusion 



In this paper the four-loop verification of identity (j3Up has been made for a special 
class of Feynman diagrams. The proof, presented in Ref. [H], is not applicable for 
these diagrams, because structure of the diagrams sum is more complicated and there 
are additional restrictions to possible positions of the gauge field external lines. Due 
to the confirmation of the considered identity, obtained here, this identity seems to 
be also valid for arbitrary diagrams. Its origin seems to be a symmetry of the theory. 
However, this symmetry surely is not the gauge invariance. 
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To conclude we stress again, that the especial importance of the obtained identity 
is that it essentially restricts the form of divergences in the considered theory. The 
matter is that it is possible to reduce the expression for the two-point Green function 
of the gauge field to the expression, depending only on the two-point Green function of 
the matter superfield and LHS of this identity, using only the gauge invariance. Thus 
the considered identity removes all undefined expressions from the two-point Green 
function of the gauge field and relates it with the two-point Green function for the 
matter superfield. 



Acknowledgements. 

This work was partially supported by Russian Foundation for Basic Research (grant 
No 05-01-00541). 



A Explicit expressions for the functions G, / and F 



Below we present expressions for the functions G, f and F. The lower index shows 
in what order of the perturbation theory (for diagrams of the considered class) the 
corresponding function was calculated. 
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fi(q 2 ) = 0; 



(42) 



Hg 2 



1 4 f d 4 kd 4 l 
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(An expression for F 3 is not presented, because it is too large.) 



(45) 
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